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A search for neutron-antineutron (n — n) oscillation, a process with \AB\ — 2 that has been predicted by 
right-left (R-L) symmetric gauge theories, was undertaken using the 24.5 x lO"^"^ neutron-yrs exposure of Super- 
Kamiokande I, in an analysis that included the significant sources of experimental uncertainties. No evidence 
for n — n oscillation was found, the lower limit of the lifetime for neutrons bound in ^^0 was determined to be 
1.89 X 10^^ yrs at the 90% confidence level (C.L.), and the corresponding limit for the oscillation time for free 
neutrons was calculated to be 2.44 x 10* sec using a theoretical suppression factor of 1.0 x 10^'' sec^^. 

PACS numbers: 



Searches for Baryon Number Violation (BNV) are moti- 
vated by various Grand Unification Theories (GUTs) and by 
the observation first put forward by A. Sakharov in 1967 that 
BNV, C, CP violation, and non-equilibrium thermodynamics 
are needed to explain the Baryon Number Asymmetry of the 
Universe (BAU) 1 1]. In the three decades since the first large 
(5 kton) water Cherenkov and fine grained nucleon decay ex- 
periments appeared in the early 1980's, the impetus for study- 
ing \B — L\ violating interactions such as neutron-antineutron 
(n — n) oscillation and double beta decay has grown steadily. 
The possibility of n — n oscillation was first discussed by 
V. Kuzmin in 1970 and recent experimental limits on 
p e+TT^ and p — )■ DK'^ have already ruled out the simplest 
\B — L\ conserving GUTs, minimal SU(5) and the minimal 
supersymmetric version of SU(5). It has also been shown that 
the triangle anomalies involving electro-weak bosons wash 
out any \B — L\ conserving baryon asymmetry above the 10 
TeV scale. The search for \B — L\ violating reactions has be- 
come increasingly important as a potential explanation of the 
observed BAU. 

The discovery of neutrino oscillations has renewed interest 
in theories with Majorana spinors which yield B and L sym- 
metry breaking by allowing \ A{B — L)\ ^ 2 with |AL|=2, 
as in neutrino-less double beta decay and |Ai3|=2 for n — n 
oscillation. These models include a large class of supersym- 
metric and R-L symmetric SU{2)l ® SU{2)ii®SU{A)c the- 
ories. Neutron-antineutron oscillation has also been predicted 
by recent GUTs with large or small extra space-time dimen- 
sions iilllQl and by R-L symmetric theories that include the 
seesaw mechanism to generate neutrino masses ||5t]. 

An important property of 7i — n oscillation is its depen- 
dence on a six quark operator, with a mass scaling of m^^ 
instead of m^'^ as in the charged X and Y bosons that mediate 
nucleon decay in minimal SU(5). Therefore, the observation 
of a significant n — n oscillation signal at Super-Kamiokande 
would imply new physics at a scale of approximately 100 TeV, 
a factor of 10 larger than the maximum energies scheduled for 
study at the Large Hadron Collider (LHC). 

The previous best 90 % C.L. lifetime limits for bound neu- 
trons are from 1MB with 1.7 and 2.4 xlO"^^ yrs. 1,^ in two 
different analyses, Kamiokande with 4.3 x 10"^^ yrs. in oxy- 
gen [7], and Soudan 2 with 7.2 x 10^^ yrs. in iron I's]. The 
current best limit on the oscillation time of unbound neutrons 
is given by ILL (Grenoble) as 0.86 x 10^ sec H. The hfe- 
time limit for bound neutrons can be converted to the n — n 
oscillation time for free neutrons by the relationship 



Tji — n — R 



(1) 



where t„_,-j and Tn-n are the oscillation time of a free neu- 
tron and the lifetime of a bound neutron, and R is the nuclear 
suppression factor, for which theoretical estimates are given 
in the literature ||22]|23|. In this paper, we describe the re- 
sults of our neutron-antineutron oscillation search in Super- 
Kamiokande I. No positive signal was observed, but from our 
data we obtained a lower limit on the lifetime of a neutron 
in oxygen of 1.89 x 10^^ yrs (90% C.L.), which is about 4.5 
times larger than the previous highest experimental limits. 

Super-Kamiokande is a ring imaging water Cherenkov 
counter containing 50 ktons of ultra-pure water, located 
in Kamioka-town in Gifu prefecture, Japan. Descriptive 
overviews of Super-Kamiokande and technical details are 
given in the literature fill. Data were taken by Super- 
Kamiokande I from May 31st, 1996 to July 15th, 2001, using 
a fiducial volume of 22.5 kton. The n — n oscillation analysis 
described in this paper used the complete SK-I data set with a 
livetime of 1489 days. 

An antineutron is expected to annihilate quickly with one of 
the surrounding nucleons and to produce multiple secondary 
hadrons, mainly pions. The total momentum of the secon- 
daries is nearly zero, except for the Fermi momenta of the an- 
nihilating n and nucleon, and the total energy is nearly equal 
to the total mass of the two nucleons. Thus one would ex- 
pect an isotropic multi-hadron signature with a total energy of 
about 2 GeV to occur in neutron-antineutron events in Super- 
Kamiokande, prior to any nuclear interactions. 

A detailed Monte Carlo (MC) simulation was formulated 
for this study. Since the literature on (n + nucleon) annihi- 
lation in nuclei is sparse, pp and^ data from hydrogen and 
deuterium bubble chambers | T2l I Hi ifT^ were used to deter- 
mine the branching ratios for the annihilation final states. The 
branching ratios included in our simulations are displayed in 
Table H] The values of kinematic quantities in our MC were 
determined using relativistic phase-space distributions that in- 
cluded the Fermi momentum of the annihilating nucleons. Pi- 
ons and omegas produced in the n-nucleon annihilations were 
propagated through the residual nucleus. This study used a 
program originally developed for the 1MB experiment to sim- 
ulate meson-nucleon interactions The cross sections for 
the pion-residual nucleus interactions were based on an in- 
terpolation from measured pion-carbon and pion-aluminum 
cross sections to those for pion-^^O interactions. Excitation 
of the A (1232) resonance was the most important effect in 
the nuclear propagation phase. It was assumed that the pion 
and omega cross sections scaled linearly with matter density, 
a quantity that decreases as the mesons move away from the 
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annihilation point and exit from the nucleus. The Fermi mo- 
mentum of the interacting nucleon and the possibility of Pauli 
blocking were also included in our simulations. We found 
that 49% of the pions did not interact, while 24% were ab- 
sorbed and 3% interacted with a nucleon to produce an addi- 
tional pion or occasionally two more pions, and the rest of the 
pion interactions involved scattering. These interaction prob- 
abilities yielded total and charged pion multiplicities of 3.5 
and 2.2, respectively, with an average charged pion momen- 
tum of 310 MeV/c and an RMS width of 190 MeV/c. There 
was a probability of 0.56% that an event included an lcP that 
emerged from the nucleus without having decayed. The final 
states of the nuclear fragments were calculated using an al- 
gorithm based on a simulation from ORNL Iil5il . Fragments 
of the residual nucleus that contained two or more nucleons 
were not simulated in water, since most of the nucleons in 
these fragments had momenta that were below the threshold 
for inelastic hadron interactions. Therefore only free n and p 
fragments were propagated through water. 

Since the largest source of background events for n — n 
oscillation is atmospheric neutrinos, we prepared a large MC 
event sample that corresponded to an atmospheric neutrino ex- 
posure of 500 yrs. A detailed description of this simulation 
is given in reference (idn . The propagation of particles and 
Cherenkov light in the detector was modeled by a program 
that is based on the GEANT-3 [17] package. The detector ge- 
ometry, the generation and propagation of Cherenkov radia- 
tion from charged particles, and the response of the PMTs and 
data acquisition electronics were also included in our simula- 
tions. Hadron interactions were simulated using the CALOR 
package OlSll for nucleons and charged pions with p^^ >500 
MeV/c, and through a custom-made program [19] for charged 
pions of <500 MeV/c. 
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TABLE I: The branching ratios for the n+nucleon annihilations in 
our simulations. These factors were derived from pp and pd bubble 
chamber data| 12]| 13][14]. 

The trigger threshold for our search corresponds to a 5.7 
MeV electron signal and yields a trigger rate of about 10 Hz. 
The background triggers are mostly due to inherent radioac- 
tivity in the detector and to cosmic ray muon interactions. 
Most of the cosmic ray muon events are rejected by a require- 
ment that there be no significant activity in the outer detector 
of Super-Kamiokande. We also require that the reconstructed 



visible energy be greater than 30 MeV to remove the remain- 
ing low energy radioactivity. This constraint is also necessary 
to maintain the effectiveness of the event reconstruction. The 
reduction algorithms are identical to those used for the atmo- 
spheric neutrino analyses and nucleon decay searches U6.] . 

Events remaining after the reduction procedure are pro- 
cessed by the full reconstruction program, which yields an 
overall event vertex, the number of Cherenkov rings, and a 
direction, particle identification determination, and a momen- 
tum for each ring For n — n oscillation events the ver- 
tex resolution is 26 cm. Each Cherenkov ring is identified as 
being either "showering" (e,7) or "non-showering" (/i,7r,p) 
based upon its hit pattern and opening angle. The momen- 
tum is subsequently determined using the assigned particle 
type and the number of collected photo-electrons inside a cone 
with an opening half-angle of 70 degrees after corrections for 
geometric effects and light attenuation are made. Our choice 
of cone angle completely covers the Cherenkov cone in wa- 
ter, which has an opening half-angle of about 42 degrees. In 
multiple-Cherenkov ring events, we estimate and separate a 
sample of photo-electrons for each ring using an expected 
Cherenkov light distribution. 

To isolate n — n candidates we apply additional criteria to 
the reduced fully contained (EC) event sample which is used 
for our atmospheric neutrino analyses il6ll : 
(a) The number of Cherenkov rings > 1, (b) 700 MeV < 
Visible energy < 1300 MeV, (c) MeV/c < Total momen- 
tum < 450 MeV/c and (d) 750 MeV/c^ < Invariant mass < 
1800 MeV/c^. The total momentum is defined as Ptot — 

^all — rings — )■ 



I where is the reconstructed momentum 

vector of the i-th ring. The invariant mass is defined to be 

E all — rings I 2 
i yVi 



while the total energy is defined as 
where to,; 



Etot ~ 2_ji ~ V Pi + n^i 7 where to,; is the mass of i- 
th ring assuming that showering and non-showering rings are 
from 7 rays and charged pions, respectively. The selection cri- 
teria given above were optimized to maximize the ratio e/^/b, 
where e is the signal detection efficiency and b is the number 
of background events. 

Distributions of the four reconstructed kinematic variables 
are displayed in Figure [1] Our atmospheric neutrino MC in- 
cludes the effects of to Vr neutrino oscillations, with mix- 
ing parameters of (sin^ 26, Ato2)=(1.0, 2.1 x lO"^ eV^) ob- 
tained from our previous publication on this topic [16]. Fig- 
ure [T] (b) shows the visible energy distribution after an event 
selection criterion based on the number of Cherenkov rings 
was imposed. The average visible energy for n~n MC oscil- 
lation events is about 700 MeV, much smaller than twice the 
nucleon mass. This decrease is mainly due to loses in hadronic 
energy as the result of scattering and absorption interactions 
with nuclear matter in the water 

Application of the event selection criteria (a)-(d) yielded 24 
candidate events, a detection efficiency of 12.1% from the re- 
maining n — n MC events in the final sample, and 24. 1 back- 
ground events as estimated for the 1489 days of Super-K I. 
Final event samples are displayed as scattered dots inside of 
the box-shaped scatterplots that apply criteria (c) and (d) in 
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FIG. 1 : Distributions of the kinematic variables subject to additional 
criteria at each of the reduction steps described in the text including: 
(a) Number of rings, (b) Visible energy, (c) Total momentum, and (d) 
Invariant mass. Circles indicate data points with statistical error-bars, 
solid black lines and dashed lines represent the atmospheric neutrino 
MC and the n — n MC, respectively. Light vertical lines with arrows 
indicate the event selection criteria. 



Figure. 12] Of the remaining background: 57.8% comes from 
multi-hadron production by Deep Inelastic Scattering (DIS), 
31.4% from single pion production via resonances, and 5.5% 
from single 77 and K meson production via resonances. The 
remaining 5.2% of the background results from Charged Cur- 
rent (CC), Quasi-Elastic (QE) scattering and Neutral Current 
(NC) elastic scattering accompanied by energetic knocked-out 
nucleons in the water. 

Uncertainties in the detection efficiency, exposure, and 
background rates are given in Table The major sources 
of errors in the detection efficiency are the results of uncer- 
tainties in the models for propagation of pions and omega 
mesons through the residual nucleus. In particular, the error 
due to uncertainties in the 7r-nucleon cross section is 20.0%, 
as estimated from the tt-^^O scattering data shown in [6] 
and by comparing results from two independent nuclear in- 
teraction programs originally developed by 1MB [6] and by 
Kamiokande 101 (6.1%). The uncertainty from the (n -1- nu- 
cleon) annihilation branching ratios is estimated to be 4.6% 
by comparing different MC results based on variations in 
the assumed branching ratios (25]|2^. The 0.6% asymme- 
try in the detector gain and the 2.0% difference in the en- 
ergy scale between the data and the MC affect the momen- 
tum cut and contribute 1 .7% and 0.4% uncertainties, respec- 
tively. By comparing real data with MC events we estimate 
that the uncertainty in the Cherenkov ring finding procedure 
is 2.2%. The detection efficiency is estimated to have a to- 
tal uncertainty of 22.9% while the exposure uncertainty was 
found to be 3.0% through an estimate of the error in the cal- 



FIG. 2: Total momentum vs. the invariant mass after applying se- 
lection criteria (a) and (b) on the PC sample: (a) Signal MC, (b) 
atmospheric neutrino MC, and (c) Data. 



culation of the effective fiducial volume. The uncertainties in 
the atmospheric neutrino flux, the atmospheric neutrino en- 
ergy spectrum, and the components of the atmospheric neu- 
trino flux contribute 7.8% to the systematic uncertainty of the 
background rate. The DIS cross section uncertainty in the low 

region yields a large contribution to the uncertainty in the 
background rate, which is estimated to be 14.1% by compar- 
ing the results from two different parameterizations of the par- 
ton distribution functions ll20ll ll2lll . The uncertainty in the en- 
ergy scale, the lack of perfect uniformity in the PMT gain, 
and the uncertainty in the Cherenkov ring finding contribute a 
total of 14.9% to the systematic errors. The total systematic 
uncertainty in the background rate is estimated to be 23.7%. 

No significant excess was found in our full 1489 day Super- 
Kamiokande I data set. The lower limit on the lifetime of 
a neutron bound inside an oxygen nucleus due to n — n os- 
cillation was calculated from the 24 observed candidate and 
24.1 expected background events. All of the systematic un- 
certainties were included in the limit calculation by employing 
Bayesian statistical methods 12711 as follows: 



P{T\n) = A 



xP{T)P{X)P{e)P{b)dXdedb. 



(2) 



The above normalization constant A in Eq.© was determined 
by imposing the constraint P(r\n)dr — 1 where F, 
A = NT, and e are the true values of the event rate, expo- 
sure, and detection efficiency, respectively, for n — n oscilla- 
tion, and b is the true mean number of background events. 
P(F),P(A),P(e),P(6) in Eq.© are the prior probability 
density functions, which we assume are Gaussian distribu- 
tions for A, e, and b. P{T) is a flat distribution for F > 
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Detection Efficiency 


Sources 


Uncertainty (%) 


Fermi momentum of nucleons 


6.2 


Annihilation branching ratio of n+nucleons 


4.6 


n propagation modeling 


6.1 


TT-nucleon cross section in the nucleus 


20.0 


Energy scale 


1.7 


Asymmetry of detector gain 


0.4 


Cherenkov ring finding 


2.2 


Total 


22.9 


Exposure 


Sources 


Uncertainty (%) 


Fiducial volume 


3.0 


Detector livetime 


<0.1 


Total 


<3.0 


Background Rate 


Sources 


Uncertainty(%) 


{Ue + i?e)/(fM + J^m)' / Tatio 


0.1, 1.0 


Up/down, horizontal/vertical flux ratios 


<1 


K/n ratio 


3.1 


Neutrino energy spectrum (<lGeV, >lGeV) 


6.1,3.6 


Neutrino cross sections 




QE 


5.8 


l-TT production 


2.4 


DIS 


14.1 


coherent tt productions 


<1 


CC/NC cross section ratio 


5.0 


Axial vector mass in QE and In production 


0.6 


Fermi momentum for QE 


<1 


TT propagation in ^''0 


4.1 


Energy scale 


4.8 


Asymmetry of detector gain 


0.5 


Cherenkov ring finding 


14.1 


Total 


23.7 



by Dover et al. ll23ll . the corresponding limit for a free neutron 
at the 90% C.L. is 



> 2.44 X 10** sec 



(6) 



compared to r„_fi > 0.86 x 10® sec as was measured by the 
ILL/Grenoble experiment |@] . 

The lifetime and deduced oscillation time limits for bound 
and free neutrons from previous n — n oscillation searches us- 
ing water Cherenkov detectors and iron calorimeters are com- 
pared with the new Super-Kamiokande I results of Eqs. (|5]) 
and (|6]l in Table Hill These earlier experiments did not include 
systematic uncertainties in the computation of their final re- 
sults, as was done for Super-Kamiokande I. 

The new results are the first in over twenty years from a 
large water Cherenkov counter, an improvement over the pre- 
vious best water Cherenkov limits for bound neutrons (from 
Kamiokande) by a factor of more than 4.5. Our results may 
also help to constrain theories of grand unification, in partic- 
ular those involving R-L symmetry, the seesaw mechanism, 
and extra dimensions. 
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(10''^neutron-yr) 


245.5 


21.9 


5.0 


3.0 


3.2 


Efficiency(%) 


12.1 


18.0 


30.0 


33.0 


50.0 


Candidates 


24 


5 








3 


B.C. 


24.1 


4.5 


2.5(2.1) 


0.9 




T^-n (10^'yr) 


1.89 


0.72 


0.65 


0.43 


0.24 


Suppression factor 












(10^^sec~i) 


1.0 


1.4 


1.4 


1.0 


1.0 


T„-fi(10'*sec) 


2.44 


1.3 


1.2 


1.2 


0.88 



TABLE II: Systematic uncertainties in efficiency, exposure, and 
background rate. 

and for negative F. Finally, the limit of the neutron lifetime 
in oxygen, Tn-n, with the inclusion of systematic uncertain- 
ties for n — n oscillation at the 90% C.L., is determined from 
eqs.OJ, and (01) as follows: 

Tn — n ^ l^lirait^ (3) 

C.L. = ( P[T\n)dT (4) 

^0 

where C.L. =0.9. The calculated result is 



TABLE 111: A comparison of the Super-Kamiokande 1 results with 
those of previous n — n experiments using bound neutrons 
llQI . The equivalent oscillation times for free neutrons were deduced 
from those obtained directly from bound neutrons using the nuclear 
suppression factor R and the expression Tn^n = R ■ r^_fi, as in 
equation (1). 
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T„_fi > 1.89 X 10^^ yrs. (5) 

This lifetime limit for a bound nucleon is converted to the 
n ~ h oscillation time of a free neutron using Eq. ([T]). Under 
the assumption that R in ^^O is 1.0 x lO'^'' sec^^, as predicted 
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